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Abstract 

It should be regarded that the confirmation of the maximum oscillation in neutrino oscillation through L/E 
analysis by Super-Kamiokande is a logical consequence of their establishment on the existence of neutrino 
oscillation through the analysis of the zenith angle distribution for atmospheric neutrino events. In the 
present paper (Part 1) with the computer numerical experiment, we examine the assumption made by 
Super-Kamiokande Collaboration that the direction of the incident neutrino is approximately the same as 
that of the produced lepton, which is the cornerstone in their L/E analysis, and we find this approximation 
does not hold even approximately. In a subsequent paper (Part 2), we apply the results from Figures 16, 17, 
18 and 19 to L/E analysis and conclude that one cannot obtain the maximum oscillation in L/E analysis 
in the single ring muon events due to quasi-elastic scattering reported by Super-Kamiokande which shows 
strongly the oscillation pattern from the neutrino oscillation. 
PACS: 13.15.-l-g, 14.60.-z 

Keywords: Super-Kamiokande Experiment, QEL, Numerical Computer Experiment, Neutrino Oscillation, 
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1. Introduction 

In principle, it is pretty difficult to specify the 
oscillation parameters in neutrino oscillation reli- 
ably from the cosmic ray experiments (atmospheric 
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neutrinos), even if they really exist, because in the 
cosmic ray experiment by its nature, one cannot de- 
termine the direction of the incident neutrino which 
plays a decisive role in the analysis of neutrino os- 
cillation and may be determined in the accelerator 
experiments even if the experimental errors are so 
large due to their geometries. 

On the other hand, according to the results ob- 
tained from the Super-Kamiokande Experiments on 
atmospheric neutrinos, it is said that oscillation 
phenomena have been found between two neutrinos, 
and v T . Published reports on the confirmation 
to the oscillation between the neutrinos, and v T , 
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and the history foregoing these experiments will be 
critically reviewed and details are in the following: 

(1) During 1980's Kamiokande and 1MB observed 
the smaller atmospheric neutrino flux ratio 
v^/ve than the predicted value 

(2) Kamiokande found anomaly in the zenith angle 
distribution Q. 

(3) Super-Kamiokande found v^-v r oscillation Q 
and Soudan2 and MACRO confirmed the 
Super-Kamiokande result [J]. 

(4) K2K, the first accelerator-based long baseline 
experiment, confirmed atmospheric neutrino 
oscillation [5j. 

(5) MINOS's precision measurement gives the con- 
sistent results with Super-Kamiokande ones@. 

It is well known that Super-Kamiokande Collab- 
oration examined all possible types of the neutrino 
events, such as, say, Sub-GeV e-like, Multi-GeV e- 
like, Sub-GeV //-like, Multi-GeV /i-like, Multi-ring 
Sub-GeV ^-like, Multi-ring Multi-GeV ^-like, Up- 
ward Stopping Muon Events and Upward Through 
Going Muon Events. In other words, all possi- 
ble interactions by neutrinos, such as, elastic and 
quasielastic scatterings, single-meson production 
and deep scattering are considered in their analy- 
ses. As the results of them, all topologically dif- 
ferent types of neutrino events lead to the uni- 
fied numerical oscillation parameters, say, Am 2 = 
2.4 x l(T 3 eV 2 and sin 2 26 = 1.0 Q. 

Really, these parameters are obtained from the 
unified analysis of the zenith angle distributions of 
various neutrino events. AS for reliability of the en- 
ergy estimation which plays a decisively important 
role in the survival probability of a given flavor, the 
qualities of the events range in various grades from 
the coarse to the fine. Generally, it is impossible 
to estimate the individual energy of the muon due 
to the neutrino interaction among Partially Con- 
tained Events due to their stochastic characters. In 
the multi-ring events, energy estimations has much 
uncertainty, even if they belong to Fully Contained 
Events. Moreover, it is absolutely impossible to es- 
timate energy of the individual event in Upward 
Going Muon Events. 

Consequently, in spite of such a fatal defect for 
the detection of the neutrino oscillation, which is 
inherent in the cosmic ray experiments as men- 
tioned above, if one goes ahead and challenges to 



prove the existence of the neutrino oscillation in 
cosmic ray experiments, then, one should focus on 
analysis of the most reliable events with the high- 
est quality exclusively from the experimental point 
of view, discarding all other physical events which 
have uncertainties more or less compared with the 
events with the highest quality. Such the most reli- 
able events among all candidates to be analyzed for 
Super-Kamiokande Collaboration are Fully Con- 
tained Events among the single ring muon events 
due to the quasi elastic scattering (QEL). These 
events occupy the majority among Fully Contained 
events observed by Super-Kamiokande Collabora- 
tion. By the definition of Fully Contained Events 
of the single ring events due to QEL, one can dis- 
criminate kinds of neutrinos and can estimate their 
transferred energies because of their confinement in 
the detector as well as their zenith angles which 
are strongly connected with their emitted angles 
in the neutrino interactions concerned. If one can 
find some clear indication on the neutrino oscilla- 
tion from the analysis of the most reliable events 
with the highest quality, then, one can surely ex- 
pect any corroboration on the neutrino oscillation 
from the analysis of neutrino events with poorer 
qualities, if they really exist. This is our motiva- 
tion for the performance of our computer numerical 
experiment. 

Here, it should be emphasized strongly that 
Super-Kamiokande Collaboration put a fundamen- 
tal assumption through their whole analyses. This 
assumption, however, is never self-evident and, 
therefore, it should be carefully examined. This as- 
sumption is that the directions of the incident neu- 
trinos are approximately the same as those of emit- 
ted leptons, which we abbreviate as the SK assump- 
tion on the direction in the present paper. This as- 
sumption should be recognized as the cornerstone 
for their neutrino oscillation analysis throughout all 
their works, which links with the survival probabil- 
ity of a given flavor, P{v^ — > v^) given by 



1 - sinHe ■ sin 2 : (1.27 'Am 2 'L v /E v ), (1) 



where sin 2 29 = 1.0 and Am 2 = 2.4 x lO~ 3 eV 2 0- 
The most important result among the achieve- 
ments on neutrino oscillation made by Super- 
Kamiokande Collaboration is the finding of the 
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maximum oscillation in neutrino oscillation, be- 
cause it is the ultimate result in the sense that they 
observe the oscillation pattern itself directly in neu- 
trino oscillation. 

Therefore, it is desirable to carry out careful L/E 
analysis for the single ring muon events due to QEL 
(the events with the highest quality) among Fully 
Contained Events exclusively in order to obtain a 
definite conclusion around the neutrino oscillation, 
focusing on the validity (or invalidity) of the SK 
assumption on the direction. Consequently, let us 
examine the SK assumption on the direction. This 
is the main theme of the present paper (Part 1). In 
a subsequent paper (Part 2), we discuss the applica- 
tion of the result obtained in Part 1 to the analysis 
for the single ring muon events due to QEL among 
Fully Contained Events, concluding the invalidity 
of the SK assumption on the direction through L/E 
analysis. 

In order to avoid any misunderstanding toward 
the SK assumption on the direction, we reproduce 
this assumption from the original SK papers and 
their related papers in italic. 

[A] Kajita and Totsuka [§] state B 

"However, the direction of the neutrino 
must be estimated from the reconstructed 
direction of the products of the neutrino 
interaction. In water Cheren-kov detec- 
tors, the direction of an observed lepton 
is assumed to be the direction of the neu- 
trino. Fig. 11 and Fig. 12 show the esti- 
mated correlation angle between neutrinos 
and leptons as a function of lepton mo- 
mentum. At energies below 400 MeV/c, 
the lepton direction has little correlation 
with the neutrino direction. The correla- 
tion angle becomes smaller with increasing 
lepton momentum. Therefore, the zenith 
angle dependence of the flux as a conse- 
quence of neutrino oscillation is largely 
washed out below 400 MeV/c lepton mo- 
mentum. With increasing momentum, the 
effect can be seen more clearly." 

[B] Ishitsuka @ statefi 

" 8.4 Reconstruction of L v 

Flight length of neutrino is determined 
from the neutrino incident zenith angle, 



2 see page 101 of the paper concerned. 
3 see page 138 of the paper concerned. 



although the energy and the flavor are 
also involved. First, the direction of 
neutrino is estimated for each sample 
by a different way. Then, the neu- 
trino flight lenght is calclulated from the 
zenith angle of the reconstructed direction. 

8.4-1 Reconstruction of Neutrino Di- 
rection 

FC Single-ring Sample 

The direction of neutrino for FC single- 
ring sample is simply assumed to be the 
same as the reconstructed direction of 
muon. Zenith angle of neutrino is recon- 
structed as follows: 

cos Ql ec = cos 6 M (8.17) 

,where cos0^ ec and cosO M are cosine of 
the reconstructed zenith angle of neutrino 
and muon, respectively." 

[C] Jung, Kajita et al. [l(| state 0: 

"At neutrino energies of more than a 
few hundred MeV, the direction of the 
reconstructed lepton approximately repre- 
sents the direction of the original neu- 
trino. Hence, for detectors near the sur- 
face of the Earth, the neutrino flight dis- 
tance is a function of the zenith direction 
of the lepton. Any effects, such as neu- 
trino oscillations, that are a function of 
the neutrino flight distance will be man- 
ifest in the lepton zenith angle distribu- 
tions." 

As clarified from Figures 11, 12 and 26 in Q and 
Figure 8.12 in [9}, Super-Kamiokande Collaboration 
know well the existence of non-negligible scattering 
angles in the neutrino interactions. Nevertheless, 
the reason why they put the fundamental postulate 
on the direction mentioned above is surmised to be 
due to their recognition that statistically enough 
accumulation of the events leads to the cosf?,, = 
cos 8 ^ as a whole after mutual cancellation of the 
scattered angles. 

As already stated, among all possible events to be 
analyzed the most important events from which we 
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can extract neutrino oscillation parameters definitly 
are undoubtedly the single ring muon (electron) 
events which are generated in the detector and ter- 
minate in it, because these events give the essential 
information for clear interpretation on neutrino os- 
cillation, if it really exists, namely, the kinds of the 
neutrinos, the transferred energies to the charged 
particles and their directions. The single muon 
events are generated mainly from the quasi-elastic 
scattering(QEL). If the existence of neutrino oscil- 
lation is verified definitely in the analysis of single 
ring muon events among Fully Contained Events 
under the SK assumption on the direction, then 
again we can say, we expect the corroborations even 
in the analysis of other types of the interactions 
with poorer qualities. Therefore, the analysis for 
the single ring muon events among Fully Contained 
Events is decisively important compared with the 
analysis of other types of the neutrino events. Con- 
sequently, first of all, let us start to examine the 
validity on the SK assumption on the direction. 

Our paper is organized as follows. In section 2, 
we treat the differential cross section for QEL in 
the stochastic manner as exactly as possible and 
obtain the zenith angle distributions of the emit- 
ted leptons, cos^, for given neutrinos with definite 
zenith angles, taking account of the azimuthal an- 
gles of the emitted leptons in QEL. As a result of it, 
we show that the SK assumption on the direction 
does not hold any more for the incident neutrinos 
with smaller energies. 

In section 3, we give the correlation between 
cos6 v (p} and cosO^ips or the correlation between 
L„( P ) and and mention the reason why the 

SK assumption on the direction does not hold. 



2. Single Ring Events among Fully Con- 
tained Events which are Produced by 
Quasi Elastic Scsattering. 

2.1. Differential cross section of quasi elastic scat- 
tering and spreads of the scattering angles 

As stated in Introduction, the finding of observa- 
tion of the maximum oscillation in the L/E analysis 
is the ultimate verification of the finding of neu- 
trino oscillation by Super-Kamiokande. For the ex- 
amination of the Super-Kamiokande's assertion, we 
analyze L/E distribution of the single ring events 
among Fully Contained Events. 

In order to examine the validity of the SK as- 
sumption on the direction, we consider the single 



ring events due to the following quasi elastic scat- 
tering(QEL): 
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The differential cross section for QEL is given as 
follows ll| . 
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The signs + and — refer to ^V e ) and Cw e ) for 
charged current (c.c.) interactions, respectively. 
The Q 2 denotes the four momentum transfer be- 
tween the incident neutrino and the charged lepton. 
Details of other symbols are given in 111. 

The relation among Q 2 , E u rp\, energy of the in- 
cident neutrino, Ei, energy of the emitted charged 
lepton (muon or electron or their anti-particles) and 
8 S , scattering angle of the emitted lepton, is given 
as 

Q 2 = 2E v{p) E t (l-cosd a ). (4) 
Also, energy of the emitted lepton is given by 
Q 2 



Ei — E v (p) - 



2M 



(5) 



Now, let us examine magnitude of the scattering 
angle of the emitted lepton in a quantitative way, 
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Figure 1: Relation between the energy of the muon and its 
scattering angle for different incident muon neutrino ener- 
gies, 0.5, 1, 2, 5, 10 and 100 GeV. 



as this plays a decisive role in determining the ac- 
curacy of direction of the incident neutrino, which 
is directly related to reliability of the zenith angle 
distribution of single ring muon (electron) events 
in the Super-Kamiokande Experiment. By using 
Eqs. to ([5]), we obtain the distribution func- 
tion for scattering angle of the emitted leptons and 
the related quantities by a Monte Carlo method. 
The procedure for determining scattering angle for 
a given energy of the incident neutrino is described 
in Appendix A. Figure Q] shows this relation for 
muon, from which we can easily understand that 
the scattering angle 9 S of the emitted lepton (muon 
here) cannot be neglected. For a quantitative exam- 
ination of the scattering angle, we construct the dis- 
tribution function for 8 S of the emitted lepton from 
Eqs. ([3]) to ([5]) by using the Monte Carlo method. 

Figure [5] gives the distribution function for 9 S of 
the muons produced in the muon neutrino interac- 
tions. It can be seen that the muons produced from 
lower energy neutrinos are scattered over wider an- 
gles and that a considerable part of them are scat- 
tered even in backward directions. Similar results 
are obtained for anti-muon neutrinos, electron neu- 
trinos and anti-electron neutrinos. 

Also, in a similar manner, we obtain not only 
the distribution function for the scattering angle of 
the charged leptons, but also their average values 
< 8 S > and their standard deviations a s . Table 1 
shows them for muon neutrinos, anti-muon neutri- 
nos, electron neutrinos and anti-electron neutrinos. 
From Table 1, it seems to be clear that the scatter- 
ing angles could not be neglected, taking account of 
the fact that the frequency of the neutrino events 
with smaller energies is far larger than that of the 
neutrino events with larger energies due to high 
steep of the neutrino energy spectrum. However, 
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Figure 2: Distribution functions for the scattering angle of 
the muon for muon-neutrino with incident energies, 0.5 , 
1.0 and 2 GeV. Each curve is obtained by the Monte Carlo 
method (one million sampling per each curve). 
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Figure 3: Schematic view of the zenith angles of the charged 
muons for different zenith angles of the incident neutrinos, 
focusing on their azimuthal angles. 
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Super-Kamiokande Collaboration assume that the 
direction of the neutrino is approximately the same 
as that of the emitted lepton even for the neutrino 
events with smaller energies, as cited in the three 
passages mentioned above 

B, S], Hoi- However, 
this assumption has never been verified by Super- 
Kamiokande Collaboration. 



2.2. Influence of the azimuthal angle in QEL over 
the zenith angle of single ring events 

In the present subsection, we examine the effect 
of the azimuthal angles, <p, of emitted leptons over 
their own zenith angles, for given zenith an- 

gles of the incident neutrinos, Qu(v)) hi QEL, which 
was not be considered in the detector simulation 
carried by Super-Kamiokande Collaboration^. The 
influence of this effect over the zenith angle cannot 
be neglected particularly in horizontal-like neutrino 
events. 

For three typical cases (vertical, horizontal and 
diagonal), Figure [3] gives a schematic representa- 
tion of the relationship between, #„( P ), the zenith 
angle of the incident neutrino, and (9 Sl 0), a pair 
of scattering angle of the emitted muon and its az- 
imuthal angle. Zenith angle of the emitted muon is 
derived from 9 V ^ and (9 S , (f>) by (A. 6) as shown in 
Appendix A. 

From Figure |3]-a, it can been seen that the zenith 
angle of the emitted lepton is not influenced 

by its (f> in the vertical incidence of the neutrinos 
= 0°), as it must be. From Figure |3J-b, how- 
ever, it is obvious that the influence of (f> of the 
emitted leptons on their own zenith angle is the 
strongest in the case of horizontal incidence of the 
neutrino {0 V ^ = 90°). Namely, one half of the 
emitted leptons are recognized as upward going, 
while the other half is classified as downward going 
ones. The diagonal case ( Q v (p) — 43°) is intermedi- 
ate between the vertical and the horizontal. In the 
following, we examine the cases for vertical, hori- 
zontal and diagonal incidence of the neutrinos with 
different energies, say, E v ^ = 0.5 GeV, E v (p\ = 1 
GeV and 2S„(p) = 5 GeV, as the typical cases. 



J Throughout this paper, we measure the zenith angles 
of the emitted leptons from the upward vertical direction of 
the incident neutrino. Consequently, notice that the sign of 
our direction is opposite to that of the Super-Kamiokande 
Experiment ( our cos0„( S ) = - cos0„( S ) in SK) 



2.3. Dependence of the spread of the zenith angle 
for emitted leptons on the energy of emitted 
leptons for different incident directions of the 
neutrinos with different energies 

The detailed procedure for our Monte Carlo sim- 
ulation is described in Appendix A. We give the 
scatter plots between the fractional energy of the 
emitted muons and their zenith angle for a definite 
zenith angle of the incident neutrino with different 
energies in Figures S] to |B1 In Figure @] we give the 
scatter plots for vertically incident neutrinos with 
different energies 0.5, 1 and 5 GeV. In this case, 
the relations between the emitted energies of the 
muons and their zenith angles are unique, which 
comes from the definition of the zenith angle of the 
emitted lepton. However, the densities (frequencies 
of the event number) along each curves are differ- 
ent in position to position and depend on the ener- 
gies of the incident neutrinos. Generally speaking, 
densities along the curves become greater toward 
cos# M (^) = 1. In this case, cos is never influ- 
enced by the azimuthal angle in the scattering by 
the definition^]. 

On the contrast, it is shown in Figure [5] that 
the horizontally incident neutrinos give the widest 
zenith angle distribution for the emitted muons 
with the same fractional energies due to the effect 
of the azimuthal angles. The lower the energies of 
the incident neutrinos are, the wider the spreads of 
the scattering angles of emitted muons # M become, 
which leads to wider zenith angle distributions for 
the emitted muons. As easily understood from Fig- 
ure H3 the diagonally incident neutrinos give the in- 
termediate zenith angle distributions for the emit- 
ted muons between those for vertically incident neu- 
trinos and those for horizontally incident neutrinos. 

It should be noticed from the figures that the 
influence of the azimuthal angle in QEL over the 
cosines of the zenith angle for the incident neu- 
trino becomes effective in the more inclined neu- 
trino, even if the scattering angle due to QEL is not 
so large. The effect of the azimutal angle in QEL is 
not taken into account in the Monte Carlo simula- 
tion by Super-Kamiokande Collaboration, because 
their Monte Carlo simulation is a detector simula- 
tion and this effect cannot be taken into account by 
their nature. 



6 The zenith angles of the particles concerned are mea- 
sured from the vertical upward direction. 
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Figure 4: The scatter plots between the fractional energies of the produced muons and their zenith angles for vertically incident 
muon neutrinos with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Figure 5: The scatter plots between the fractional energies of the produced muons and their zenith angles for horizontally 
incident muon neutrinos with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 1000 for each case. 



(a) 



(b) 



(c) 



1 

0.8 
0.6 
0.4 
0.2 



-0.2 
-0.4 
-0.6 
-0.8 



muon 










cos9 v =C 
— E^=€ 


.731 (9 V = 
).5GeV 


=43") 














■., • J . 









































































0.2 



0.4 0.6 

E„ / E v 



0.8 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 



muon 














































• ''• 






































■•' cose 


Y =0.731 


(9 V =43°) 






E 


v =1Ge\ 





0.2 0.4 0.6 0.8 



1 

0.8 
0.6 
0.4 
0.2 



-0.2 
-0.4 
-0.6 
-0.8 



muon 
















• '...•> 


.■ •z&r. 


































































cose 


v =0.731 


(9 V =43°) 






E 


v =5Ge\ 


1 



0.2 



0.4 

E„ 



0.6 0.8 



/E„ 



Figure 6: The scatter plots between the fractional energies of the produced muons and their zenith angles for diagonally 
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Figure 7: Zenith angle distribution of the muon for the vertically incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stands for the corresponding ones under the SK assumption on 
the direction. 
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Figure 8: Zenith angle distribution of the muon for the horizontally incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stands for the corresponding ones under the SK assumption on 
the direction. 



(a) 



(b) 



mu 

coi 


Dn r 


eut 
0.7 


rino 

3K<jv=4 


3") 










a\ 


E v = 
/era 


0.5 

ge= 


0.1 


39 






s 


K 






s 


d.= 


0.5£ 


56 










































































.. . ■ 



















-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS9„ 



6 
5 

# 4 
co 
o 
o 

5 3 
z 

a 

2 
1 





mu 

co: 


Dn r 

i9 v = 


eut 
0.7 


'ino 
31(6 


v=4 


3°) 










a\ 


E v = 
/era 


1Gt 

ge= 


3V 

0.4 


32 






S 


K 






s 


d.= 


0.46 


53 









































































-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS9„ 



(c) 



muon reut 
co!;9 v =0.7Sl(9 v =4t3°) 



=5Gi5V 
a^erage=0.7|5 



d.=3.103 



SK 



-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS9„ 



Figure 9: Zenith angle distribution of the muon for the diagonally incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stands for the corresponding ones under the SK assumption on 
the direction. 



Table 1: The average values < 9 B > for scattering angle of the emitted charged leptons and their standard deviations a s for 
various primary neutrino energies E v tj}\. 



E v(p) (GeV) 


angle 
(degree) 










0.2 


<o s > 


89.86 


67.29 


89.74 


67.47 




Os 


38.63 


36.39 


38.65 


36.45 


0.5 


<0 S > 


72.17 


50.71 


72.12 


50.78 




Os 


37.08 


32.79 


37.08 


32.82 


1 


<0 S > 


48.44 


36.00 


48.42 


36.01 




Os 


32.07 


27.05 


32.06 


27.05 


2 


<0s> 


25.84 


20.20 


25.84 


20.20 




Os 


21.40 


17.04 


21.40 


17.04 


5 


<0s> 


8.84 


7.87 


8.84 


7.87 




Os 


8.01 


7.33 


8.01 


7.33 


10 


<0s> 


4.14 


3.82 


4.14 


3.82 




Os 


3.71 


3.22 


3.71 


3.22 


100 


<0s> 


0.38 


0.39 


0.38 


0.39 




Os 


0.23 


0.24 


0.23 


0.24 



2.4- Zenith angle distribution of the emitted leptons 
for different incident directions of the neutri- 
nos with different energies 

In Figures [7] to GO we express Figures 2] to [5] in a 
different way. We sum up muon events with differ- 
ent emitted energies for given zenith angles. As the 
result of it, we obtain frequency distribution of the 
neutrino events as a function of cos0^ for different 
incident directions and different incident energies of 
neutrinos. 

In Figures [TJa) to[7fc), we give the zenith angle 
distributions of the emitted muons for the case of 
vertically incident neutrinos with different energies, 
say, E v = 0.5, 1 and 5 GeV. 

Comparing the case for 0.5 GeV with that for 
5 GeV, we understand the big contrast between 
them as for the zenith angle distribution. The scat- 
tering angle of the emitted muon for 5 GeV neu- 
trino is relatively small (See, Table 1), so that the 
emitted muons keep roughly the same direction as 
their original neutrinos. In this case, the effect of 
their azimuthal angle on the zenith angle is also 
smaller. However, in the case for 0.5 GeV which is 
the dominant energy for single ring muon events in 
the Super-Kamiokande, there is even a possibility 
for the emitted muon to be emitted in the back- 
ward direction due to the larger angle scattering, 
the effect of which is enhanced by their azimuthal 



angle. 

The most frequent occurrence in the backward 
direction of the emitted muon appears in the hor- 
izontally incident neutrino as shown in Figs. 8(a) 
to 8(c). In this case, the zenith angle distribution 
of the emitted muon should be symmetrical with 
regard to the horizontal direction. Comparing the 
case for 5 GeV with those both for ~0.5 GeV and 
for ~1 GeV, even 1 GeV incident neutrinos lose al- 
most the original sense of the incidence if we mea- 
sure it by the zenith angle of the emitted muon. 
Figures HJJa) to[Sfc) for the diagonally incident neu- 
trinos tell us that the situation for diagonal case lies 
between the case for the vertically incident neutri- 
nos and that for horizontally incident ones. SK in 
the figures denotes the SK assumption on the direc- 
tion of incident neutrinos. From the Figures [7| a) to 
M,c), it is clear that the scattering angles of emitted 
muons influence their zenith angles, which is en- 
hanced by their azimuthal angles, particularly for 
more inclined directions of the incident neutrinos. 

2.5. Zenith Angle Distribution of Fully Contained 
Events for a Given Zenith Angle of the Inci- 
dent Neutrino, Taking Their Energy Spectrum 
into Account 

In the previous sections, we discuss the relation 
between the zenith angle distribution of the inci- 
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Figure 10: Zenith angle distribution of fi~ and fi + for v and v for the incident neutrinos with the vertical, horizontal and 
diagonal directions, respectively (see Figure 3). The overall neutrino spectra at Kamioka site are taken into account. The 
sampling number is 10000 for each case. SK stand for the corresponding ones under the SK assumption on the direction. 



dent neutrino with a definite energy and that of the 
emited muons produced by the neutrino for the dif- 
ferent incident direction. In order to apply our in- 
spection around the uncertainty of the SK assump- 
tion on the direction for Fully Contained Events to 
the real situation, we must consider the effect of 
the energy spectrum of the incident neutrino. The 
Monte Carlo simulation procedures for this purpose 
are given in Appendix B. 

In Fig. [TUJ we give the zenith angle distributions 
of the sum of and /i _ for a given zenith angle 

of Vjx and z/^, taking into account different primary 
neutrino energy spectra for different directions at 
Kamioka site. SK in the figures denotes the SK as- 
sumption on the direction. From the figures, it is 
clear that the SK assumption on the direction does 
not hold. Namely, we can conclude that the scat- 
tering angle of the emitted muons acompanied by 
their azimuthal angles influence their zenith angle 
distribution for all directions. 

3. Super-Kamiokande Assumption on the 
Direction and the Real Relation between 
cos9 v (p} and cosO^p) and the correspond- 
ing relation between L v (p^ and 

3.1. Correlation between cos 6 ' u andcosO^ 

Now, we extend the results for the definite zenith 

angle obtained in the previous sections to the case 

in which we consider the zenith angle distribution 

of the incident neutrinos totally. 

Here, we examine the real correlation between 

cos6 v and cos0 M , by performing the exact Monte 

Carlo simulation. 



The detail for the simulation procedure is given 
in Appendix C. 

In order to obtain the zenith angle distribution of 
the emitted leptons for that of the incident neutri- 
nos, we divide the range of cosine of the zenith angle 
of the incident neutrino into twenty regular inter- 
vals from cos V — to cos0„ = 1. For the given in- 
terval of cos Q Vl we carry out the exact Monte Carlo 
simulation, and obtain the cosine of the zenith an- 
gle of the emitted leptons. 

Thus, for each interval of cos#„, we obtain the 
corresponding zenith angle distribution of the emit- 
ted leptons. Then, we sum up these corresponding 
ones over all zenith angles of the incident neutri- 
nos and we finally obtain the relation between the 
zenith angle distribution for the incident neutrinos 
and that for the emitted leptons. 

In a similar manner, we could obtain between 
cos 9p and cos 6p for anti-neutrinos. The situation 
for anti-neutrinos is essentially the same as that for 
neutrinos. 

Super-Kamiokande Collaboration have made the 
comprehensive analysis of the all possible types of 
the neutrino events on the neutrino oscillations 0] 
for 1489.2 live days. We restrict our analysis to the 
neutrino events due to QEL to make uncertainties 
around the interpretation small as possible. 

In Figure 111! we give the correlation diagram be- 
tween cos^j, and cos(9 p for 1489.2 live days with- 
out oscillation. Here, we classify the correlations 
with regard to different incident neutrino ener- 
gies. It is clear from Figure [TT] that the neu- 
trino events for E v > 5 GeV exist on the line 
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Figure 11: Correlation diagram between cos0„ and cos 9^ for null oscillation for different neutrino energy regions. 
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Figure 12: Correlation diagram between cos0„ and cos 9^ for oscillation for different neutrino energy regions. 
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Figure 13: Correlation diagram between cos and cos M for the events which exist in null oscillation but disappear due to 
oscillation for different neutrino energy regions. 



cos 9 V = cos 9^ roughly, and the neutrino events for 
2 < E v < 5 GeV exist along the line cos 9 V — cos M 
with explicit deviations. Namely, in higher neutrino 
energy regions, the SK assumption on the direction 
roughly holds. However, if we pay our attention 
to lower energies, for 1 < E v < 2 GeV, we easily 
understand that their distribution deviates largely 
from the line cos 9 V — cos# M , showing clearly that 
the SK assumption on the direction does not hold. 

Consequently, in sum, we can conclude that the 
SK assumption on the direction does not hold as a 
whole, taking account of the fact that lower energy 
neutrino events occupy the majority in the detector 
due to the strong steepness of the incident neutrino 
energy spectra. In Figure [T^l we give the correla- 
tion diagram between cos^ and cos 6^ for 1489.2 
live days with oscillation. The Monte Carlo simu- 
lation with oscillation is carried out with the use 
of the rejection method with regard to the survival 
probability for a given flavor, which is based on 
the results obtained by the Monte Carlo simulation 
without oscillation 0. 

The regions in the correlation diagrams in Fig- 



7 We can carry out the Monte Carlo simulation for the 
case with oscillation quite independently from that without 
oscillation, of course. Howerver, the adoption of the rejection 



ures [TT1 and [T2l are divided into the four parts, each 
of which have their own physical meaning. The 
neutrino events in the first sector where cos 9 V > 
and cos 9^ > consist of the upward neutrinos and 
upward muons. Namely, this denotes such a situa- 
tion that the upward neutrinos emit the muons in 
forward directions. The neutrino events in the sec- 
ond sector where cos 9 V < and cos 9^ > consist 
of the downward neutrinos and the upward muons. 
Namely, this denotes such a situation that down- 
ward neutrinos emit the muons in backward direc- 
tions. The neutrino events in the third sector where 
cos 9 V < and cos 9^ < consist of the downward 
neutrino and the downward muons. Namely, this 
denotes such a situation that the downward neutri- 
nos emit the muons in forward directions. The neu- 
trino events in the fourth sector where cosf?,, > 
and cos 9^ < consist of the upward neutrinos and 
the downward muons. Namely, this denotes such a 
situation that the upward neutrinos emit the muons 
in backward directions. 

Here, the downward muons are generated in three 
different manners. In the first case, they are gen- 
erated by downward neutrinos. In the second case, 



method make clearer the effect of the survival probability for 
the given flavor. 
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they are generated by upward neutrinos, owing to 
the backscattering in QEL. And in the third case, 
they are generated by either upward or downward 
horizontal-like neutrinos. In the last case, either 
upward or downward horizontal-like neutrinos may 
produce downward muons due to the accidental az- 
imuthal effect in QEL. (see, Figure 3b). The upward 
muons are generated in three different manners sim- 
ilarly by both upward and downward neutrinos. 

It is clear from Figure [11] that events in the first 
sector are symmetrical to those in the third sector, 
while events in the second sector are symmetrical 
to those in the forth sector. As we expect the same 
incident neutrino fluxes in both downward and up- 
ward directions in the case of no oscillation, these 
symmetries are should be hold, as it must be. On 
the other hand, it is clear from Figure IT2"1 that sym- 
metries found in Figure [TT] are lost any more. This 
comes from that the upward neutrino flux is smaller 
than that of the downward neutrino in the case of 
the existence of oscillation. 

In Figure 1131 we give the events which exist in 
the case of no oscillation but disappear due to os- 
cillation obtained through the Monte Carlo oper- 
ational procedure by the survival probability of a 
given flavor ( Eq.(l) ). It is simply obtained by 
the subtraction of the events in Figure [T^] from the 
events in Figure [TT] 

We can obtain several interesting results from 
Figure [TJ|] The first is that we do not find dis- 
appeared (rejected) events in the downward neu- 
trino events due to oscillation. This reason is 
very simple and clear. It is due to the extreme 
low probability of oscillation for downward neu- 
trino with the statistics considered on the occa- 
sion of the selection of the specified neutrino os- 
cillation paremeters by Supcr-Kamiokande Collab- 
oration, say Am 2 = 2.4 x 10~ 3 eV 2 . When we adopt 
this value one order of magnitude larger than SK's, 
Am 2 = 2.4 x 10~ 2 eV 2 , we surely expect the events 
concerned for the downward neutrinos. Further- 
more, when we adopt this value one order of magni- 
tude smaller than SK's, Am 2 = 2.4 x 10~ 4 eV 2 , then 
we surely expect different distribution of the events 
concerned even for the upward neutrinos, keeping 
the absence of the events concerned for downward 
neutrinos. The second is that almost disappeared 
(rejected) events belong to the first sector (upward 
neutrinos — > upward muons). The third is that 
we cannot neglect disappeared (rejected) neutrino 
events which belong to the forth sector (upward 
neutrino — > downward muon). 




Figure 14: Schematic view of relations among L„, L M , 8 S 
and <h a . 
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Figure 15: The procedure for our numerical experiment for 
obtaining L M from a given L v . 
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3.2. Correlation between L„ and 

In the previous subsection, we verify that the SK 
assumption on the direction does not hold from the 
examination on the correlatin between cos 6 V and 
cos# M . 

This SK assumption on cosf?^ and cos 9^ is log- 
ically equivalent to the statement that L v is ap- 
proximately the same as in L/E analysis, where 
L v denotes the distance on the incident neutrino 
from the interaction point of the neutrino events to 
the intersection of the Earth surface toward its ar- 
riving direction and denotes the corresponding 
distance on the emitted muon. Consequently, if our 
indication on the invalidity of the SK assumption 
on the direction on cos9 u and cos 8^ is correct, the 
same conclusion should be drawn from the correla- 
tion between L„ and L^. In the present subsection, 
we directly examine the validity of the implicit SK 
assumption that L v is approximated by L M , taking 
consideration of the neutrino energy spectrum at 
the Super-Kamiokande site. The relation between 
L v and is given in Figure [14] Also, we show 
the procedure to determine for a given L v in 
Figure [15] 

L v and are functions of the direction cosine 
of the incident neutrino, cosO v , and that of emitted 
muon, cosO^, respectively and they are given as, 



Lv-Lu correlation for FC 



L v = R g x(rsKCos6„+ 



SK 1 



?9v 



1 - r% K ) 



R g x(rsKCOs9^+^Jr 



2 SK cos 2 e il 



i 



SKI 



(6-1) 
(6-2) 



where R g is the radius of the Earth and r$K — 
1 — Dsx/Rg, with the depth, Dsk, of the Super- 
Kamiokande Experiment detector from the surface 
of the Earth. It should be noticed that the L„ and 
are regulated by both the energy spectrum of 
the incident neutrino and the production spectrum 
of the muon (QEL in the present case). Conse- 
quently, their mutual relation is influenced by either 
the absence of oscillation or the presence of oscilla- 
tion which depend on the combination of oscillation 
parameters. 

In Figures Ql)] and [XT] we give the correlation di- 
agram between L v and for single ring muon 
events among Fully Contained Events for the 1489.2 
live days in both the absence and the present of 
neutrino oscillation which corresponds to the actual 
Super-Kamiokande Experiment [7j. In Figure llSl we 
give the subtraction of Figure [T7] from Figure [TH] 
which corresponds to Figure 1131 exactly. In Fig- 
ures [16] to [TSl blue points denote neutrino events 
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Figure 16: Correlation diagram for L v and without oscil- 
lation for 1489.2 live days. The blue points and orange points 
denote neutrino events and ani-neutrino events, respectively. 
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Figure 17: Correlation diagram for L u and L M with oscilla- 
tion for 1489.2 live days. The blue points and orange points 
denote neutrino events and ani-neutrino events, respectively. 
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while orange points denote anti-neutrino events. 
Throughout all correlation diagrams in the present 
paper, blue points and orange ones have the same 
meaning as shown in these Figures. The aggregates 
of the (anti-) neutrino events which correspond to 
definite combinations of L v and L M are essentially 
classified into four groups in the following, and they 
correspond to four sections with regards to cosr?^ 
and cos 9^ as shown in Figures [TT] to [13J 

Group A is defined as the aggregate for neutrino 
events in which both L v and are rather small. It 
denotes that the downward neutrinos produce the 
downward muons with smaller scattering angles. In 
this case, energies of the produced muons are near 
those of the incident neutrinos due to smaller scat- 
tering angles. 

Group B is defined as the aggregate for neutrino 
events in which both L v and are rather large. 
It denotes that the upward neutrinos produce up- 
ward muons with smaller scattering angles. In this 
case, the energy relation between incident neutri- 
nos and produced muons is essentially the same as 
in Group A in Figure [T51 because the flux of the 
upward neutrino events is symmetrical to that of 
downward neutrino events in the absence of neu- 
trino oscillation. 

Group C is defined as the aggregate for neutrino 
events in which L v are rather small and are 
rather large. It denotes that the downward neu- 
trinos produce the upward muons by the possible 
effect resulting from both backscattering and az- 
imuthal angle in QEL. In this case, energies of the 
produced muons are smaller than those of the inci- 
dent neutrinos due to larger scattering angles. 

Group D is defined as the aggregate for the neu- 
trino events in which L u are rather large and are 
rather small. It denotes that the upward neutrinos 
produce the downward muons. The energy relation 
between the incident neutrinos and the produced 
muons is essentially the same as in Group C in the 
absence of neutrino oscillation (Figure IT61) . 

Summarized from the mentioned above, we can 
say that there exist the symmetries between Group 
A and Group B, and also between Group C and 
Group D, which reflect the symmetry between the 
upward neutrino flux and the downward neutrino 
one for null oscillation. 

In Figure [T71 we give the correlation between L v 
and under their neutrino oscillation parameters, 
say, Am 2 = 2.4 x l(T 3 eV 2 and sin 2 28 = 1.0 Q. 
In the presence of neutrino oscillation, the prop- 
erty of the symmetry which holds in the absence 



of neutrino oscillation (see (Group A and Group 
B) and/or (Group C and Group D) in Figure ITS)) 
is lost due to different incident neutrino fluxes in 
the upward direction and the downward one. If we 
compare Group A with Group B, the event num- 
ber in Group B (upward v — > upward /i) is smaller 
than that in group A (downward v — > downward 
jti), which comes from smaller flux of the upward 
neutrinos. The similar relation between Group C 
(downward v — » upward \x ) and Group D (upward 
v — > downward fi) holds in Figure 1171 

Summarizing the characteristics among Groups 
A to D in the Figures [T|)] and [T71 we can conclude 
that (Group A and Group B) and (Group C and 
Group D) are in symmetrical situations in the ab- 
sence of neutrino oscillation, while such a symmet- 
rical situation is lost in the presence of neutrino 
oscillation. Also, it is clear from Figures [H)] and 
1171 that L v pa Lfj,, namely the SK assumption on 
the direction, does not hold both in the absence of 
neutrino oscillation and in the presence of neutrino 
oscillation even if statistically. 

Here, it should be noticed that the approxima- 
tion of L v w does not hold completely in the 
regions C and D. The event numbers in Group C 
and Group D could not be neglected among the 
total event number concerned. In these regions, 
neutrino events consist of those with backscatter- 
ing and/or neutrino events in which the neutrino 
directions are horizontally downward (upward), but 
their produced muons turn upward (downward) re- 
sulting from the effect of azimuthal angles in QEL. 

In the L/E analysis made by Super-Kamiokande 
Collaboration, the reconstruction of the direction 
of the incident neutrino from the direction of the 
emitted muon is very simple, say, 

cos Q r u ec = cos 0^ 

and there is nothing more0. 

Of course, as we pointed out already, Super- 
Kamiokande Collaboration know the existence of 
the non-negligible scattering angles of the neutrino 
interaction. However, they don't consider it in 
their analysis. The concept of the "reconstruction" 
by Super-Kamiokande Collaboration denote one-to- 
one correspondence between the direction of the in- 
cident neutrino and the direction of the emitted lep- 
ton. However, the correlations between them shown 



8 following "cos©£ ec = cos © M (8.17)" , Ishitsuka states 
"cos©£ ec and cos © M are cosines of the reconstructed zenith 
angles of neutrino and muon, respectively." (see text in page 
3) . 
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Figure 19: The correlation diagram between E v and Eu (a) without oscillation, (b) with oscillation and (c) for the events 
which exist in null oscillation but disappear due to oscillation, for 1489.2 live days, respectively. The solid line denotes the 
polynomial expression by Super-Kamiokande Collabolation. 



in Figures [16] and [17] make it impossible to recon- 
struct the direction of the incident neutrino from 
the direction of emitted muon, even if one takes the 
scattering angle into account, because these corre- 
lations are the reflection of the stochastic charac- 
ters in the processes concerned which deny unique 
relation between them. Again, we should empha- 
size and confirm that "reconstruction of the direc- 
tion of the incident neutrino" in L/E analysis of 
Super-Kamiokande Collaboration is just equal to 
the adoption of the direction of the emitted muon 
due to the neutrino interaction. 

3.3. The correlation between E v and 

Here, we examine the transformation of E v from 
Eft adopted by Super-Kamiokande Collaboration. 
The validity problem is closely related to the SK 
assumption on the direction through the survival 
probability for a given flavor. 

Super-Kamiokande Collaboration estimate E v 
from Efj,, the visible energy of the muon, from 
their Monte Carlo simulation, by the following 
equation^ (see page 135 of the paper concerned) : 

Eu,sk ^E^x(a + bxx + cxx 2 + dx x 3 ), (7) 

where x = log^E^). 

The idea that E v can be approximated as the 
polynomial of E^ means that there should be the 



unique relation between E v and E^. However, in 
the light of stochastic characters inherent in both 
the incident neutrino energy spectrum and the pro- 
duction spectrum of the muon, such a treatment 
is not suitable theoretically, which may kill a real 
correlation effect between the incident neutrino en- 
ergy and the emitted muon energy. In Figures Q1J1 
we give the correlation between E v and E^ to- 
gether with that obtained from the polynomial ex- 
pression by Super-Kamiokande Collaboration un- 
der their neutrino oscillation parameters and their 
incident neutrino energy spectrum [12]. It is clear 
from the figure that the part of the lower energy 
incident neutrino deviates largely from the approxi- 
mated formula, which reflects explicitly the stochas- 
tic character of QEL. 

In Figure [19] (b), we give the correlation be- 
tween E v and E^ with oscillation. Glancing at Fig- 
ures [in] (a) and (b), we cannot recognize the dif- 
ference between them, because too much events are 
marked on the figure for their discrimination. Then 
in Figure IT91 (c), we give the subtraction between 
them, as shown in Figures [13] and [18] It is easily 
understood that there is the visible difference be- 
tween them. Furthermore, it is easily understood 
from Figures Q1J] that the polynomial expression is 
not suitable for the description on the mutual rela- 
tion between E v and E^ , standing on the stochastic 



16 



point of view. 

Finally, it is necessary to mention to the qual- 
itative difference between the correlation between 
L v and L M , and the correlation between E v and 
Efj,, from the degree of their influence on the SK 
assumption on the direction. As clarified in Fig- 
ures [11] and [T7J it is impossible to approximate as 
L v nearly equal to [L v « L^). On the other 
hand, the approximation of E v by Eq.(7) does not 
introduce fatal error finally, although its treatment 
is not theoretically suitable in the sense of the lack 
of stochastic characters on the physical processes. 
Namely, what influences on the L/E analysis essen- 
tially is the correlation between L u and L M , but not 
the correlation between E v and E^. 

4. Conclusion 

Since one cannot measure L v and E u , one is 
forced to utilize L M and E^ in the L/E analysis 
in place of them. Then, Super-Kamiokande Col- 
laboration assume that the direction of the inci- 
dent neutrino is approximately the same as that of 
the emitted lepton, namely the SK assumption on 
the direction and E v can be estimated from some 
polynomial formula of the variable E^ in their L/E 
analysis. However, it is clear from Figures H2] and 
IT3l and/or Figures H"6l and [171 that the SK assump- 
tion on the direction does not hold even approxi- 
mately and the transformation of E^ into E v is not 
uniquely. 

Consequently, we can conclude that L v cannot 
be reconstructed by i M , but E v can be done if we 
allow considerable uncertainties. 

In the subsequent paper (Part 2), we apply the 
results from Figures [T1J] to QjJ] to L/E analysis 
and conclude that one cannot obtain the maxi- 
mum oscillation in L/E analysis reported by Super- 
Kamiokande which shows strongly the oscillation 
pattern due to neutrino oscillation. 
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Appendix A. Monte Carlo Procedure for 
the Decision of Emitted Ener- 
gies of the Leptons and Their 
Direction Cosines 

Here, we give the Monte Carlo Simulation pro- 
cedure for obtaining the energy and its direction 
cosines, (l r ,m r ,n r ), of the emitted lepton in QEL 
for a given energy and its direction cosines, (I, m, n), 
of the incident neutrino. 

The relation among Q 2 , E u ^\, the energy of the 
incident neutrino, E^j\, the energy of the emitted 
lepton (muon or electron or their anti-particles) and 
6 S , the scattering angle of the emitted lepton, is 
given as 

Q 2 = 2E u(p) E m (l ~ cos# s ). (AT) 
Also, the energy of the emitted lepton is given by 
Q 2 



2M 



(A2) 



Procedure 1 

We decide Q 2 from the probability function for 
the differential cross section with a given -E„( P ) 
(Eq. (J3]) in the text) by using the uniform random 
number, £, between (0,1) in the following 



Q 2 

£= / Pi { l)(E„ {P) ,Q 2 )dQ 2 , (A3) 
where 

&vl(l){E v {v), Q 2 



dQ 2 




dQ 2 . 
(AT) 



From Eq. (AT), we obtain Q 2 in histograms to- 
gether with the corresponding theoretical curve 
shown in Figure IATI The agreement between the 
sampling data and the theoretical curve is excellent, 
which shows the validity of the utlized procedure in 
Eq. (A-3). 

Procedure 2 

We obtain E t ^ from Eq. (A-2) for the given £v( P ) 
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Q 2 /Q 2 n 



Figure A.l: The reappearance of the probability function for 
QEL cross section. Histograms are sampling results, while 
the curves concerned are theoretical ones for given incident 
energies. 




(/ , m , n ) 



( m, n) \ 



and Q 2 thus decided in Procedure 1. 
Procedure 3 

We obtain cos6* s , cosine of the the scattering angle 
of the emitted lepton, for E e ^ thus decided in the 
Procedure 2 from Eq. (A-l) . 

Procedure 4 

We decide <f>, the azimuthal angle of the scattered 
lepton, which is obtained from 



2tt£. 



(A-5) 



Here, £ is a uniform random number between (0, 
I)- 

As explained schematically in the text (see Figure [3] 
in the text), we must take account of the effect 
due to the azimuthal angle <f> in QEL to obtain the 
zenith angle distribution both for Fully Contained 
Events and Partially Contained Events correctly. 

Procedure 5 

The relation between direction cosines of the in- 
cident neutrino, (^(p) , Tn u t^ , n v i p \ ) , and those of 
the corresponding emitted lepton, (£ r ,m r ,n r ), for a 
certain 9 S and <f> is given as 




VP 



V -VPT> 



VP 



o 



sin# s cos(j 
sinf? s sin</ 
cos# s 



where n v ( B ) 



(A-6) 

= cos9 v ^) 1 and n? = cosQi. Here, 6i is 



the zenith angle of the emitted lepton. 

The Monte Carlo procedure for the determi- 
nation of 0£ of the emitted lepton for the parent 
(anti-)neutrino with given 6 u rp) and E v ip\ involves 
the following steps: 



Figure A. 2: The relation between the direction cosine of the 
incident neutrino and that of the emitted charged lepton. 



We obtain (£ r ,m r7 n r ) by using Eq. (|A-6|) . The 
n r is the cosine of the zenith angle of the emitted 
lepton which should be contrasted with ro„, that of 
the incident neutrino. 

Repeating the procedures 1 to 5 just mentioned 
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above, we obtain the zenith angle distribution of 
the emitted leptons for a given zenth angle of the 
incident neutrino with a definite energy. 



In the SK analysis, instead of Eq. (IA-6|1 . they 
assume n r — n v <p\ uniquely for E^-\ > 400 MeV. 

Appendix B. Monte Carlo Procedure to 
Obtain the Zenith Angle of 
the Emitted Lepton for a 
Given Zentith Angle of the 
Incident Neutrino 

The present simulation procedure for a given 
zenith angle of the incident neutrino starts from 
the atmospheric neutrino spectrum at the opposite 
site of the Earth to the SK detector. We define, 

,no— osc 

(£^( P ), i, cos#„(p)), the interaction neu- 
trino spectrum at the depth t from the SK detector 
for the case no oscillation in the following way, 

^Vint,no_osc(^iy(p), t, COS0„(p)) = 

dt \ 



1 - 



M{Eu{p),tx,px) 



x ••■ x 1 



dt 



j in i Pn) 



(B-l) 



Here, N ap (E v f 9 \, cos is the atmospheric 

(anti-) neutrino spectrum for the zenith angle at the 
opposite surface of the Earth. And Xi(E u rp\,ti, pi) 
denotes the mean free path from QEL for the neu- 
trino (anti neutrino) with the energy E v ip\ at the 
distance, U, from the opposite surface of the Earth, 
where pi is there density. 

In the presence of oscillation, neutrino energy 
spectrum correponding to (B-l) is given as, 

-^intjOsc [Ejj^i/^ , t, cos9 v ^)) 

= A r int,no_osc(^(p),C0s6' Iy (p)) X P(v^ -> V^) 

(B-2) 

Here, Piy^ — > v^) is the survival probability of a 
given flavor, such as u^, and it is given by 

P(yn -> v v) = 

1 - sin 2 26 ■ sin 2 '(1.27 'Am 2 ] L V /E U ), (B-3) 



where sin 2 28 = 1.0 and Am 2 = 2.4 x 10~ 3 eV 2 
obtained from Super-Kamiokande Collaboration^. 

The procedures of the Monte Carlo Simulation 
for the incident neutrino(anti neutrino) with a 
given energy, E v ^ , whose incident direction is 
expressde by (l,m,ri) are as follows. 

Procedure A 

For the given zenith angle of the incident neutrino, 
@u(p), w e formulate, N plo (E u ^,t, cos 6 ' v (u))dEv(p), 
the production function for the neutrino flux to pro- 
duce leptons at the Kamioka site in the following 

= a t{E) {E v (p))Ni n t (E V ( S ) , t, cos0 J/ (y) )dB 1/ ( P ) , 

(B-4) 



where 



, fQmax da t(l) (E v(p) ,Q 2 ) 
<?ew\E v (p)) - ^2 dQ . 



(B-5) 



Each differential cross section above is given by Eq. 
© in the text. Here, we simply denote the inter- 
action energy spectrum as N- m t(E v rp},t,cos0 v (p\) , 
irrespective of the absence or the presence of oscil- 
lation. 

Utilizing, £, the uniform random number between 
(0,1), we determine E v r p \, the energy of the in- 
cident neutrino in the following sampling procedure 



Pd(E v (p),t,cos v (p)(f))dE v ^), 



(B-6) 



where 



Pd(Ev(p) ! COS Q V (y~) 

_ N pro (E v ^,t, cos 9 v (p) )dE v ^) 

-(C), max 

N pro {E u (p) , t, cos 9 v (p) )dE v (p) 

(B-7) 

In our Monte Carlo procedure, the reproduction 
of, Pd{E u ( 9 -\,t, cos6 v /p\)dE v tp\, the normalized dif- 
ferential neutrino interaction probability function, 
is confirmed in the same way as in Eq. (A- 4). 
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Procedure B 

For the (anti-)neutrino concerned with the energy 
of E v rp\, we sample Q 2 utlizing £, the uniform 
random number between (0,1). The Procedure B 
is exactly the same as Procedure 1 in Appendix A. 

Procedure C 

We decide, S , the scattering angle of the emitted 
lepton for given £/„(p) and Q 2 . Procedure C is ex- 
actly the same as in the combination of Procedures 
2 and 3 in Appendix A. 

Procedure D 

We randomly sample the azimuthal angle of the 
charged lepton concerned. The Procedure D is 
exactly the same as in Procedure 4 in Appendix A. 

Procedure E 

We decide the direction cosine of the charged 
lepton concerned. Procedure E is exactly the same 
as Procedure 5 in Appendix A. 

We repeat Procedures A to E until we reach the 
desired trial number. 



P„(cos 6 v (p) 



I{cosO v{9) ) 



/ J(cos0„(p))d(cos0„(p)) 
Jo 



(C-3) 



Procedure B 

For the sampled d(cos#„( P )) in Procedure A, we 



sample E v i&\ from Ea. (IC-4[) by using £, the uniform 
randum number between (0,1) 



Ppro {E v {9) , COS l/ (p))dS„(p) , (C-4) 



where 



Ppro{E u (v), t, COS u (p) )dE v ( 9 ) — 

N pro (E L ,^,t, cos S !/ (p))dB I/ (p) 
Np IO (E u (p} , t, cos $i,(p) 

(C5) 



Appendix C. Correlation between the 
Zenith Angles of the Incident 
Neutrinos and Those of the 
Emitted Leptons 

Procedure A 

By using, N pro ( E v ^ , t, cos 9 U ^ )dE v ^ , which is 
defined in Eq. (|B-4[) . we define the spectrum for 
cos v [y\ in the following. 

i"(cos0„(p))d(cos0„( P )) = 
d(cos^(p)) x 

x / N pTO (E v ^, t, cos^( P ))d£' l/ (p). 

(C-1) 

By using Eq. (IC-2[) and £, a sampled uniform ran- 
dom number between (0,1), then we can determine 
cos Q v {p) from the following equation 

/•COS 

£= / P„(cos^ (p) )d(cos^ ( p ) ), (C-2) 

Jo 

where 



Procedure C 

For the sampled E u fo\ in Procedure B, we sample 
Efj,(p,) from Eqs. (|A-2[) and (|A-3j) . For the sampled 
-E„(p) and E^/p,), we determine cos# s , the scatter- 
ing angle of the muon uniquely from Eq. (|A- 1|) . 

Procedure D 

We determine, <fi, the azimuthal angle of the 
scattering lepton from Eq. (| A-5|) by using £, an 
uniform randum number between (0,1). 

Procedure E 

We obtain cosfl^p) from Eq. (|A-6j) . As the result, 
we obtain a pair of (cos^rp), cos 9^)) through 
Procedures A to E. Repeating the Procedures A 
to E, we finally obtain the correlation between the 
zenith angle of the incident neutrino and that of the 
emitted muon. 
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